Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense.

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories.

The genetic factors underlying growth traits differ over time points or stages. However, most current studies of phenotypes at single time points do not capture all loci or explain the genetic differences underlying growth trajectories. Hybrid Liriodendron exhibits obvious heterosis and is widely cultivated, although its complex genetic mechanism underlying growth traits remains unknown. A genome-wide association study (GWAS) is an effective method for elucidating the genetic architecture by identifying genetic loci underlying complex quantitative traits. In the present study, using a GWAS, we identified robust loci associated with growth trajectories in hybrid Liriodendron populations. We selected 233 hybrid progenies derived from 25 crosses for resequencing, and measured their tree height (H) and diameter at breast height (DBH) for 11 consecutive years; 192 972 high-quality single nucleotide polymorphisms (SNPs) were obtained. The dynamics of the multiyear single-trait GWAS showed that year-specific SNPs predominated, and only five robust SNPs for DBH were identified in at least three different years. Multitrait GWAS analysis with model parameters as latent variables also revealed 62 SNPs for H and 52 for DBH associated with the growth trajectory, displaying different biomass accumulation patterns, among which four SNPs exerted pleiotropic effects. All identified SNPs also exhibited temporal variations in effect sizes and inheritance patterns potentially related to different growth and developmental stages. The haplotypes resulting from these significant SNPs might pyramid favorable loci, benefitting the selection of superior genotypes. The present study provides insights into the genetic architecture of dynamic growth traits and lays a basis for future molecular-assisted breeding.

Genome-wide identification of XTH genes in Liriodendron chinense and functional characterization of LcXTH21.

Liriodendron chinense is a relic tree species of the family Magnoliaceae with multiple uses in timber production, landscape decoration, and afforestation. L. chinense often experiences drought stress in arid areas. However, the molecular basis underlying the drought response of L. chinense remains unclear. Many studies have reported that the xyloglucan endotransglucosylase/hydrolase (XTH) family plays an important role in drought stress resistance. Hereby, to explore the drought resistance mechanism of L. chinense, we identify XTH genes on a genome-wide scale in L. chinense. A total of 27 XTH genes were identified in L. chinense, and these genes were classified into three subfamilies. Drought treatment and RT-qPCR analysis revealed that six LcXTH genes significantly responded to drought stress, especially LcXTH21. Hence, we cloned the LcXTH21 gene and overexpressed it in tobacco via gene transfer to analyze its function. The roots of transgenic plants were more developed than those of wild-type plants under different polyethylene glycol (PEG) concentration, and further RT-qPCR analysis showed that LcXTH21 highly expressed in root compared to aboveground organs, indicating that LcXTH21 may play a role in drought resistance through promoting root development. The results of this study provide new insights into the roles of LcXTH genes in the drought stress response. Our findings will also aid future studies of the molecular mechanisms by which LcXTH genes contribute to the drought response.

Overexpression of the LcCUC2-like gene in Arabidopsis thaliana alters the cotyledon morphology and increases rosette leaf number.

BACKGROUND: The unique 'mandarin jacket' leaf shape is the most famous trait of Liriodendron chinense and this characteristic gives L. chinense aesthetic and landscaping value. However, the underlying regulatory mechanism of genes involved in the leaf development of L. chinense has remained unclear. METHODS: Based on transcriptome data of leaves at different developmental stages from L. chinense, we identified differentially expression genes (DEGs) functioning in leaf development. A candidate gene named LcCUC2-like (LcCUC2L) had high similarity in sequence with Arabidopsis thaliana CUC2, and used for further research. We isolated the full-length LcCUC2L gene and its promoter from L. chinense. Subsequently, we analyzed the function of the LcCUC2L gene and its promoter activity via transformation into A. thaliana. RESULTS: In this study, we found that the LcCUC2L and AtCUC2 are homologous in sequence but not homologous in function. Unlike the role of AtCUC2 in leaf serration and SAM formation, the LcCUC2L mainly regulates cotyledon development and rosette leaf number. Histochemical ß-glucuronidase (GUS) staining revealed that LcCUC2L was expressed in the cotyledons of A. thaliana seedlings, indicating that the LcCUC2L may play a role in cotyledon development. Ectopic expression of LcCUC2L resulted in long, narrow cotyledons without petioles, abnormal lamina epidermis cells and defective vascular tissue in cotyledons, and these results were consistent with the LcCUC2L expression pattern. Further analysis showed that overexpression of LcCUC2L also induced numerous rosette leaves. Also, LcCUC2L and other related genes showed a severe response in L. chinense by introducing exogenous auxin stimulation, partly revealed that LcCUC2L affects the leaf development by regulating the auxin content. CONCLUSIONS: These results suggest that LcCUC2L may play a critical role in leaf development and morphogenesis in L. chinense, and our findings provide insight into the molecular mechanisms of leaf development in L. chinense.

Transcriptomic profiling of the floral fragrance biosynthesis pathway of Liriodendron and functional characterization of the LtuDXR gene.

Floral fragrance, which has the function of attracting pollinators, is a class of volatile secondary metabolites mainly released by the secretory tissue of petals. Terpenoids are key components of floral volatile substances. Previous studies have shown that there are significant differences in the concentration and composition of volatile floral fragrances, especially terpenoids, between Liriodendron chinense and L. tulipifera. At present, the mechanism by which the synthesis of floral fragrance is regulated in Liriodendron remains unexplored. In this study, we analyzed the differentially expressed genes (DEGs) of L. chinense and L. tulipifera, and identified 130 DEGs related to terpenoid synthesis. A KEGG enrichment analysis of DEGs related to terpenoid biosynthesis revealed that the monoterpenoid biosynthesis pathway was the most significant. We cloned the LtuDXR gene from L. tulipifera using RACE technology. RT-qPCR results showed that the expression of the LtuDXR gene was the highest in the early florescence petals, indicating that the LtuDXR gene may play a role in the synthesis of volatile terpenoids. Subcellular localization showed that the LtuDXR protein is mainly localized in the chloroplast. Overexpression of LtuDXR in Arabidopsis thaliana significantly increased the plant height, DXR enzyme activity, and carotenoid content. In this study, we identified and functionally characterized LtuDXR, which is involved in terpenoid synthesis in Liriodendron. Our work lays the foundation for further exploration of the molecular mechanism by which terpenoid biosynthesis is regulated in Liriodendron.

Genetic divergence and local adaptation of Liriodendron driven by heterogeneous environments.

Ecological adaptive differentiation alters both the species diversity and intraspecific genetic diversity in forests, thus affecting the stability of forest ecosystems. Therefore, knowledge of the genetic underpinnings of the ecological adaptive differentiation of forest species is critical for effective species conservation. In this study, single-nucleotide polymorphisms (SNPs) from population transcriptomes were used to investigate the spatial distribution of genetic variation in Liriodendron to assess whether environmental variables can explain genetic divergence. We examined the contributions of environmental variables to population divergence and explored the genetic underpinnings of local adaptation using a landscape genomic approach. Niche models and statistical analyses showed significant niche divergence between L. chinense and L. tulipifera, suggesting that ecological adaptation may play a crucial role in driving interspecific divergence. We detected a new fine-scale genetic structure in L. chinense, and divergence of the six groups occurred during the late Pliocene to early Pleistocene. Redundancy analysis (RDA) revealed significant associations between genetic variation and multiple environmental variables. Environmental association analyses identified 67 environmental association loci (EALs; nonsynonymous SNPs) that underwent interspecific or intraspecific differentiation, 28 of which were associated with adaptive genes. These 28 candidate adaptive loci provide substantial evidence for local adaptation in Liriodendron. Our findings reveal ecological adaptive divergence pattern between Liriodendron species and provide novel insight into the role of heterogeneous environments in shaping genetic structure and driving local adaptation among populations, informing future L. chinense conservation efforts.

Genome-wide survey and identification of AP2/ERF genes involved in shoot and leaf development in Liriodendron chinense.

BACKGROUND: Liriodendron chinense is a distinctive ornamental tree species due to its unique leaves and tulip-like flowers. The discovery of genes involved in leaf development and morphogenesis is critical for uncovering the underlying genetic basis of these traits. Genes in the AP2/ERF family are recognized as plant-specific transcription factors that contribute to plant growth, hormone-induced development, ethylene response factors, and stress responses. RESULTS: In this study, we identified 104 putative AP2/ERF genes in the recently released L. chinense genome and transcriptome database. In addition, all 104 genes were grouped into four subfamilies, the AP2, ERF, RAV, and Soloist subfamilies. This classification was further supported by the results of gene structure and conserved motif analyses. Intriguingly, after application of a series test of cluster analysis, three AP2 genes, LcERF 94, LcERF 96, and LcERF 98, were identified as tissue-specific in buds based on the expression profiles of various tissues. These results were further validated via RT-qPCR assays and were highly consistent with the STC analysis. We further investigated the dynamic changes of immature leaves by dissecting fresh shoots into seven discontinuous periods, which were empirically identified as shoot apical meristem (SAM), leaf primordia and tender leaf developmental stages according to the anatomic structure. Subsequently, these three candidates were highly expressed in SAM and leaf primordia but rarely in tender leaves, indicating that they were mainly involved in early leaf development and morphogenesis. Moreover, these three genes displayed nuclear subcellular localizations through the transient transformation of tobacco epidermal cells. CONCLUSIONS: Overall, we identified 104 AP2/ERF family members at the genome-wide level and discerned three candidate genes that might participate in the development and morphogenesis of the leaf primordium in L. chinense.

Genome-Wide Identification and Expression Analysis of R2R3-MYB Family Genes Associated with Petal Pigment Synthesis in Liriodendron.

The MYB transcription factor family is one of the largest families in plants, and its members have various biological functions. R2R3-MYB genes are involved in the synthesis of pigments that yield petal colors. Liriodendron plants are widely cultivated as ornamental trees owing to their peculiar leaves, tulip-like flowers, and colorful petals. However, the mechanism underlying petal coloring in this species is unknown, and minimal information about MYB genes in Liriodendron is available. Herein, this study aimed to discern gene(s) involved in petal coloration in Liriodendron via genome-wide identification, HPLC, and RT-qPCR assays. In total, 204 LcMYB superfamily genes were identified in the Liriodendron chinense genome, and 85 R2R3-MYB genes were mapped onto 19 chromosomes. Chromosome 4 contained the most (10) R2R3-MYB genes, and chromosomes 14 and 16 contained the fewest (only one). MEME analysis showed that R2R3-MYB proteins in L. chinense were highly conserved and that their exon-intron structures varied. The HPLC results showed that three major carotenoids were uniformly distributed in the petals of L. chinense, while lycopene and ß-carotene were concentrated in the orange band region in the petals of Liriodendron tulipifera. Furthermore, the expression profiles via RT-qPCR assays revealed that four R2R3-MYB genes were expressed at the highest levels at the S3P/S4P stage in L. tulipifera. This result combined with the HPLC results showed that these four R2R3-MYB genes might participate in carotenoid synthesis in the petals of L. tulipifera. This work laid a cornerstone for further functional characterization of R2R3-MYB genes in Liriodendron plants.

Isolation, expression, and functional analysis of the geranylgeranyl pyrophosphate synthase (GGPPS) gene from Liriodendron tulipifera.

Terpenoids are important secondary metabolites in plants and are involved in stress responses and pollinator attraction. Geranylgeranyl pyrophosphate synthase (GGPPS) is a key synthase in the 2C-methyl-D-erythritol-4-phosphate (MEP) pathway of terpenoid synthesis, catalyzing the synthesis of diterpenoids. Liriodendron tulipifera is a nectar plant in North America. Little is known about the key genes involved in the biosynthetic pathways of terpenoids, the precursors of most compounds related to nectar, fragrance and coloring in flowers in L. tulipifera. In this study, the LtuGGPPS2 gene and its promoter (LtuGGPPS2-pro) were cloned from L. tulipifera. The results of sequence alignment showed that the LtuGGPPS2 gene is highly homologous to GGPPS genes of other plants. Subcellular localization analysis showed that the LtuGGPPS2 protein localizes to chloroplasts, suggesting that the LtuGGPPS2 gene is probably related to carotenoid and chlorophyll synthesis. Based on tissue expression profiles revealed by RT-qPCR, the expression level of the LtuGGPPS2 gene was highest in petals. These results were consistent with the changes in volatile and nonvolatile terpenoids in the flowers of L. tulipifera. GUS staining to examine the LtuGGPPS2 promoter indicated that it is responsive to hormones. Overexpression of the LtuGGPPS2 gene increased the carotenoid content and GGPPS enzyme activity in Arabidopsis thaliana, indicating that LtuGGPPS2 is the key terpenoid synthase in the flowers of L. tulipifera. Our findings lay a foundation for further functional analysis of the LtuGGPPS2 gene and deeper investigation of the terpenoid biosynthetic pathway in L. tulipifera.

Alternative Splicing Enhances the Transcriptome Complexity of Liriodendron chinense.

Alternative splicing (AS) plays pivotal roles in regulating plant growth and development, flowering, biological rhythms, signal transduction, and stress responses. However, no studies on AS have been performed in Liriodendron chinense, a deciduous tree species that has high economic and ecological value. In this study, we used multiple tools and algorithms to analyze transcriptome data derived from seven tissues via hybrid sequencing. Although only 17.56% (8,503/48,408) of genes in L. chinense were alternatively spliced, these AS genes occurred in 37,844 AS events. Among these events, intron retention was the most frequent AS event, producing 1,656 PTC-containing and 3,310 non-PTC-containing transcripts. Moreover, 183 long noncoding RNAs (lncRNAs) also underwent AS events. Furthermore, weighted gene coexpression network analysis (WGCNA) revealed that there were great differences in the activities of transcription and post-transcriptional regulation between pistils and leaves, and AS had an impact on many physiological and biochemical processes in L. chinense, such as photosynthesis, sphingolipid metabolism, fatty acid biosynthesis and metabolism. Moreover, our analysis showed that the features of genes may affect AS, as AS genes and non-AS genes had differences in the exon/intron length, transcript length, and number of exons/introns. In addition, the structure of AS genes may impact the frequencies and types of AS because AS genes with more exons or introns tended to exhibit more AS events, and shorter introns tended to be retained, whereas shorter exons tended to be skipped. Furthermore, eight AS genes were verified, and the results were consistent with our analysis. Overall, this study reveals that AS and gene interaction are mutual-on one hand, AS can affect gene expression and translation, while on the other hand, the structural characteristics of the gene can also affect AS. This work is the first to comprehensively report on AS in L. chinense, and it can provide a reference for further research on AS in L. chinense.

Cloning, Characterization and Functional Analysis of the LtuPTOX Gene, a Homologue of Arabidopsis thalianaIMMUTANS Derived from Liriodendron tulipifera.

Flower colour and colour patterns are crucial traits for ornamental species; thus, a comprehensive understanding of their genetic basis is extremely significant for plant breeders. The tulip tree (Liriodendron tulipifera Linn.) is well known for its flowers, odd leave shape and tree form. However, the genetic basis of its colour inheritance remains unknown. In this study, a putative plastid terminal oxidase gene (LtuPTOX) was identified from L. tulipifera based on multiple databases of differentially expressed genes at various developmental stages. Then, the full-length cDNA of LtuPTOX was derived from tepals and leaves using RACE (rapid amplification of cDNA ends) approaches. Furthermore, gene structure and phylogenetic analyses of PTOX as well as AOXs (alternative oxidases), another highly similar homologue in the AOX family, were used to distinguish between the two subfamilies of genes. In addition, transient transformation and qPCR methods were used to determine the subcellular localization and tissue expression pattern of the LtuPTOX gene. Moreover, the expression of LtuPTOX as well as pigment contents was investigated to illustrate the function of this gene during the formation of orange bands on petals. The results showed that the LtuPTOX gene encodes a 358-aa protein that contains a complete AOX domain (PF01786). Accordingly, the LiriodendronPTOX and AOX genes were identified as only paralogs since they were rather similar in sequence. LtuPTOX showed chloroplast localization and was expressed in coloured organs such as petals and leaves. Additionally, an increasing pattern of LtuPTOX transcripts leads to carotenoid accumulation on the orange-band during flower bud development. Taken together, our results suggest that LtuPTOX is involved in petal carotenoid metabolism and orange band formation in L. tulipifera. The identification of this potentially involved gene will lay a foundation for further uncovering the genetic basis of flower colour in L. tulipifera.